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ABSTRACT 

Here we introduce moments of visibility function and discuss how those can be used to 
estimate the power spectrum of the turbulent velocity of external spiral galaxies. We 
perform numerical simulation to confirm the credibility of this method and found that 
for galaxies with lower inclination angles it works fine. The estimator outlined here is 
unbiased and has the potential to recover the turbulent velocity spectrum completely 
from radio interferometric observations. 
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1 INTRODUCTION 

It is found that the power s pectra of the H l column density 
fluctu ations in our Galaxy (ICrovisier fc DickevlIlQsA iGreenI 
Il993h follow power law at length scales ranging a few parsec 
to few hundred parsecs. Different probes to the statistics of 
the turbulent velocity fluctuations of our Galaxy: statistics 
of centroid of velocities dEsauive l fc LazariarJl20 09l). velocit y 
coordinate spectrum or VGS ( Pogosvan fc LazariarJ 20091) . 
veloc ity channel analysis or VGA ( Lazariai^_&_^OEOS_^ad 
l200d) . spectral correlation function i Padoan et al.l l200ll ) 
show that at similar scales, the power spectrum of the 
velocity fluctuations also follow power laws. Theoretically, 
these structures are understood in term s of compressible 
fluid turbulence llElmegreen fc Scaldl200j) driven by super¬ 
novae shocks, self gravity etc. Recently, H l column density 


power spectrum f rom external dwarf (lBegum_^^_al 


Duttaet_^lJ I2OO8I) and spiral galaxies ( Dutta et al 


2ooe 


2008 


20091 . I2OI3I) are also found to follow power law at length 


scales ranging from a few parsec to ~ 10 kpc. Generat¬ 
ing mechanism of these large scal e structures are yet to 
be understood dWalker et al.l 1^141 ). Estimating the veloc¬ 
ity fluctuation power spectrum would give more clue here. 
For this purpose, VGA is the most efficient technique and 
has been used by several aut hors successfully for our g alaxy 


nas been used by several aut Uors successtully tor our galaxy 
and nearby dwarf gala xies dStanimirovic fc LazarianI 200 ll : 
iGhepurnov et all l2015l) . However, for the external spiral 
galaxies, this technique has limited use (lDuttall2015h . 

In this letter we propose to use moments of the observed 
visibilities to perform a model dependent estimate of the 
statistics of the line of sight component of turbulent veloci¬ 
ties from external galaxies. We outline an estimator for the 
power spectrum of the turbulent velocity fluctuations and 
access it’s efficacy for a noise free simulated data. Section 


(2) defines visibility moments and the estimator. Simulation 
parameters are given in section (3) and in section (4) we 
explain the procedure to use the estimator on the simulated 
data. We discuss the results and conclude in the last section 
( 5 ). 


2 POWER SPECTRUM OF VISIBILITY 
MOMENTS 

Specific intensity of radiation with rest frequency vo from 
H I gas at r = (a:, y, z) = {R, z) at a distance D and having 
te mperature T c an be written for small optical depth limit 
as (lDraindl201ll ) 


Ii9,v) = lo I dzUHlif) (j){v). 


( 1 ) 


Here v = c{vq — v)/nuo, v is the frequency of observation, 
Jo _ ^ jg jjgg q£ sight direction. R is in 

the plane of the sky, angular separation of a direction in the 
held of view with respect to the held centre is 0 — The 
line shape function (/>(v) can be written as 


(/)(v) = 00 exp 


(n - 


2cr2(T) 


( 2 ) 


were a{f) = El is the thermal velocity dispersion and 

v^(r) is the line of sight velocity component of the H l gas. 

We use radio interferometers to observe H l emissions 
from the spiral galaxies. Every antenna pair of an inter¬ 
ferometric array measures the quantity termed as visibility 
V(U, v) at a baseline U given by the projected antenna sep¬ 
aration at the plane of the sky in units of observing wave¬ 
length. Visibilities are direct Fourier transform of the specific 
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^ fci, : Boltzmann constant, mjfj : mass of hydrogen atom 
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intensity distribution, i.e, 

V{U,v) = j de (3) 

Depending on the antenna positions of the particular tele¬ 
scope array, integration time of observation and declination 
of the source, visibilities can only be measured at a finite set 
of points in the U plane. While reconstructing the image by 
performing an inverse Fourier transform of the observed vis¬ 
ibilities, above sampling results in a complica ted telescope 
beam pattern, usually known as dirty beam (iTavlor et al.l 
Il999lf . Hence, image obtained by an inverse Fourier trans¬ 
form of the visibilities is a convolution of the dirty beam 
with the sky specific intensity. Various deconvolution tech¬ 
niques are used to estimate the specific intensity from the 
image. However, all these algorithms effectively interpolate 
the visibilities in the unsampled part of the U plane. This 
gives rise to an unknown bu t correlated noise in the recon¬ 
structed image JPuttal bOllI'). Hence, any estimate of power 
spectrum of column density or line of sight velocity, evalu¬ 
ated using the reconstructed specific intensity is expected to 
have an unknown bias. 

We define the zeroth moment of the visibility function 
as follows 

Vo(t/) = J dvV{U,v). (4) 

The velocity integral above has to be performed over the 
entire spectral range of H l emission. Clearly the mo¬ 
ment Vo{U) is proportional to the the Fourier transform 
of columns density Nhi{9) ~ J dz nHi{R, z). In fact, the 
quantit y Po(i7) = U Vo(U)Vo(U)* I) has been used in lit¬ 
erature JCreenll 19931 : iDutta 8z Bharadwaill201^ ~) to estimate 
the power spectrum of the column density fluctuations of our 
Galaxy and external galaxies. We define the first moment of 
visibility function as 

Vi{U) = J dvvV{U,v) (5) 

= J de 7o J dzuHiii^Vzif). 

Here Vz(fJ, the line of sight velocity, has two components: 
v^(9), systematic rotation of the galaxy, and (9), random 
motion because of turbulence. Hence, Vz(fJ = v^{f)+vj (r). 
In case of the spiral galaxies with relatively small inclination 
angle, we may assume that the line of site component of the 
rotational velocity is independent of the distance along the 
line of sight, then 

Vi{U) = Vo(t/)®u?(i/) + xitJ), (6) 

where v^{U) is the Fourier transform of v2(r) and x{U) is 
the Fourier transform of 

xi9)=Io J dz HHiir^v'^if). (7) 


2.1 An estimator of the turbulent velocity 
fluctuations 

It is safe to assume that the systematic rotational velocity 
of the galaxy and the random turbulent velocity are not 


correlated. We define 

VxiU) = (I X0) f) = (I Vi(t7) f -V,{U)® I v^{U) 1^ ). 

( 8 ) 

The angular brackets < ... > denote ensemble averaging, 
which has to be taken over many statistical realisations of 
the observed sky. As we are provided with only one realisa¬ 
tion in reality, this can not be performed. Turbulence gen¬ 
erated fluctuations in density and velocity, particularly for 
magnetohydrodyn amic (MHD) turbulence, c an have pref- 
fered orientations dEsauivel fc LazarianI [2OO5II . However, it 
is not clear, if at the length scales of 1 kpc to 10 kpc or higher 
one would expect to see anisotropic fluctuations in neutral 
hydrogen resulting because of MHD turbulence. Here, we 
assume that we are probing homogenious and isotropic fluc¬ 
tuations only. The method presented here can be generalized 
for the anitropic case, however, observations with higher sig¬ 
nal to noise than presently available may be required. With 
this assumptions, the ensemble average in eqn. (8) can be 
replaced by averages over annular bins in U. Note that the 
quantities | Vi{U) and Vo{U) can be estimated directly 
from the observed visibilities. The line of sight component 
of the rotational velocity v^{9), inclination or position an¬ 
gle etc depends only on the galacto-centric radius. These 
quantities, and hence also | v^(U) |, can be estimated from 
the reconstructed image without any bias. It is clear from 
eqn. (7) that VxiU) contain in it informations about H l den¬ 
sity as well as H I velocity fluctuation power spectra, whereas 
the H I density fluctuation power spectra can be estimated 
independently using VoiU). 

Next we assume that the turbulent velocity and the den¬ 
sity fluctuations are not correlated and refer to the following 
two cases. 

• If we are interested in estimating the power spectra at 
length scales comparable to the scale of the galaxy’s disk, 
i.e, scales much larger than the scale height of the galaxy 
in consideration, we may approximate the galaxy as a thin 
disk. In this case, we can write 

VxiU) =Voiu)^r^iu), (9) 

where (U) = {| v^{lJ) |^) is the power spectrum of the 
turbulence velocity. It is to be noted that the quantities 
VoiU) and V^ (U) are two dimensional power spectra of the 
column density and the velocity fluctuations respectively. 
We have already discussed that VoiU) can be measured from 
the zeroth moment of the observed visibilities. Since we ex¬ 
pect all three functions in the above expression to be smooth 
and continuous, we can estimate (U) either by decon¬ 
volving VoiU) from VxiU) or by regression analysis using a 
parametric model for Vi^{U). Since we have an independent 
and exact measure of VoiU) from the zeroth moment of vis¬ 
ibilities, we can have a complete estimate of Pj (U) in this 
case. 

• On the other hand, if we are to estimate the power spec¬ 
tra at scales comparable to the scale height of the galaxy’s 
disk, we may write, 

Vx0) ^P^°,iU,O)0P^f{U,O), (10) 

where Pn^i iU, 0) and P^t iU, 0) are the three dimensional 
power spectra of nniir) and vj{r) evaluated at the zero 
line of sight wavenumber. It can be shown that Pn^j iU, 0) 
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Figure 1. Radial profile (dash line), vertical profile (dot-dash 
line) of the H I number density adopted for simulation are plot¬ 
ted against distance from the centre of the galaxy in radial and 
vertical directions. The solid black line represents the rotation 
curve used to generate v^{9). Note distances are normalised to 
the size of the simulation box. 


Figure 2. Vf (U) estimated from the simulation is shown with 
dot-dashed line. This uses the parameters given in the bottom 
left corner of the image (/3 = —2.0). A model spectrum (U) 

with 7 = —2.5 is also shown (dashed line) and it does not match 
with 'Pf{U). We also show the corresponding column density 
power spectrum (solid black line). All three spectra are nor¬ 
malised at U = 1. 


is proportional to VolU) for a thick disk. Hence, just like 
the case for the thin disk, we may use VoiU) to find the 
turbulent velocity fluctuation power spectrum from above 
relation either by deconvolution or by regression analysis. 
However, our ability to estimate the amplitude of would 
depend on the details of the morphology of the galaxy in the 
line of sight direction. 


3 SIMULATING H i OBSERVATIONS FROM A 
SPIRAL GALAXY 

A major assumption we took in deriving the expression for 
Vx{U) in the previous section is that is independent of 
the line of sight direction. In case of galaxies with lower incli¬ 
nation angle, this assumption is justified. However, for galax¬ 
ies with higher inclination angles, line of sight pass through 
different galacto-centric radius and this assumption may not 
be justified. We simulate H l observation from an external 
spiral galaxy and use it to study the effect of galaxy’s incli¬ 
nation on the turbulent velocity power spectrum estimator 
described in the previous section. Simulation o f the H l signa l 
is based on eqn. (1) and (2) and discussed in IDu tti (l2015l) 
in detail, we note the salient features here. We consider an 
exponential radial profile and a gaussian vertical profile to 
model the H l distribution in the galaxy, i.e, 

nHi{Pg) = no exp I Rg I/Ro) exp (a/zo)^^ 

[1 + fni 5nHi{rg)] , 

( 11 ) 

where Pg — {Rg, Zg) is a vector centred at the galactic centre 
with Rg at the plane of the galaxy. We assume a constant 
inclination angle i between the line of sight direction and 
Zg for simplicity. This gives a rotation mapping between the 
coordinates r and fg. Turbulent H l number density fluctu¬ 
ation 5Hi{rg) is expected to have a power law spectrum and 
we quantify it with a slope a. We assume that the tangentia l 
rotation velocity vt{Rg) follow Brandt profile jBran dliQeol) 


of the form 


Vt{Rg) = 


Rn 


Ti + 2^1 

[3 ^ 3 J 


3/2 ’ 


( 12 ) 


and use this to generate the line of sight component of the 
systematic velocity of the H l gas across the galactic disk 


..(7 - +/../W. 

Hn 


(13) 


We shall discuss the other parameters in above expressions 
shortly. 

We consider 512® cubic grids in our simulation volume 
with each cube representing an individual H l cloud. To ac¬ 
commodate the simulated galaxy’s disk inside the simula¬ 
tion volume and adopt a generic rotation curve, we consid¬ 
ered both Ro — Rv ~ 100 grid units. Assuming H l extent 
of a typical spiral galaxy is ~ 50 kpc these choice gives 
us a spatial resolution of ~ 200 pc. Simulating a thin disk 
requires relatively higher computational resources. On the 
other hand inclination angle affects both the thick and thin 
disks in similar way, we opted for simulating a thick disk 
here and adopt a value of zo — 100 grid units. We show 
these radial and vertical profiles normalised to their central 
value in Figure (1) along with the tangential velocity profile 
nor malised to up. _ 

iDutta fc Bharadwail (|2013| ) have reported that the am¬ 
plitude of the H I column dens^ fluctuation of the spiral 
galaxies in the THINGS samplqj are approximately 1/10*^ 
of the mean column density. Following this, we choose /hi = 
0.1. Line of sight velocity dispersion of H l for the galaxies 
in the THINGS sample vari es in the range ~ 5 to 20 km 
sec“® (iTamburro et al ll2009h . while the flattening value of 
the rotatio n curve are found to be in the range ~ 100 to 200 
km sec“® dde Blok et akllioOSh . Since this velocity disper¬ 
sion is mostly because of turbulence, we adopt /„ = uq/IO 
in the simulation, the actual value of vo is unimportant here. 


2 THINGS: The H I Nearby Galaxy Survey dWalter et al.ll2008h 
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Figure 3. Best fit values of /3, i.e, 0 found using the estimator to the simulated data is plotted against the input values of 0 with 
3(7 errors for four different sets of simulations. The solid line represents the ideal case oi 0 = 0. Clearly at high inclination angles the 
estimated values of /3 has relatively higher uncertainty. 


iKanekar et al.l (l201l[ l found that in our Galaxy, line of sights 
with column density > 2 x 10^'^ atoms cm^ have tempera¬ 
ture of 240 K or lower, whereas the mean column density 
of the spiral galaxies in t he TH INGS sample is > 3 x 10^° 
atoms cm^ l|Walter et al.l 120081 : IPutta fc Bharadwail 1201,311 . 
This suggest existence of large amount of cold gas in the spi¬ 
ral galaxies. In view of this, thermal velocity dispersion a is 
considered to be uo/lOO. For a uo = 200 km sec“^, this cor¬ 
responds to a gas temperature 500 K, i.e from cold neutral 
medium. 

H I column density power spectra of 18 spiral galax- 
ies in the THINGS sample are found to follow power laws 
llDutta et all l2013l l for length scales more than the scale 
height of the disk and the slope of the power law for most 
of the galaxies were found to be in the range —1.5 to —1.8. 
To see the effect of a on the efficacy of our estimator we 
perform the simulation for two different a = (—2.5,—1.5). 
We consider a nearly face on disk with i — 10° and a 
more edge on disk with i = 45°. With these choices, we 
run four sets of simulations as (a) {i = 10°, a = —1.5), 
(b) (i = 10°, a = —2.5), (c) {i = 45°, a = —1.5) and 
(d) {i = 45°,a = —2.5). Turbulent velocity component 
jr) is also taken to follow power law power sp e ctrum 
llFederrath et al.l l2009l ~l. with a slope 0. iGoldmaiJ (l2000l ~) 


pointed out that it is likely that the ISM density fluctu¬ 
ations passively follow the velocity fluctuations and hence 
they would have similar power spectra. Observations show 
that the column density fluc tuation power spectra has a 
slope of ~ —2.5 iIg reen[ jl993l l in our Galaxy and ~ —1.5 
dDutta fc Bharadwail 2013l f in case of the external galaxies. 
We simulate the H l emission for 0 varying in the range 
—0.5 to —3.5 in steps of 0.1 and generate H l specific in¬ 
tensity data cubes following eqn. (1) with a velocity reso¬ 
lution of ct/ 2. This amounts in calculating the line shape 
function for at least six points within it’s support, neces- 
sary and sufficient to avoid th e shot noise bias discussed in 
IChepurnov fc LazarianI (l2009l l. We convert these specific in¬ 
tensity maps to corresponding visibilities and use them as 
the observed visibilities V{U, v) for further analysis. In this 
letter we do not consider any measurement noise. Further, 
visibilities are estimated at grids in U space and hence the 
baseline coverage is assumed to be complete. 
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4 ESTIMATING TURBULENT VELOCITY 
FLUCTUATION FROM THE SIMULATED 
DATA 

We proceed to estimate the power spectrum of turbulent 
velocity fluctuations through the following steps. 

• We estimate the visibility moments Vo{U) and Vi(?7) 
from the visibilities V{U,v). 

• We reconstruct the image data cube from visibilities 
and estimate the rotation curve using a standard tilted ring 
fit to the data. Since, in the simulation we have kept the in¬ 
clination angle constant, we also keep it to be constant with 
radius while performing the tilted ring fit. For a real data, we 
would need to consider the uncertainties in the estimate of 
the rotation curve (and also inclination and position angles). 
At the present case these uncertainties are not considered. 
Using the rotation curve estimate we construct v^{U). 

• We bin the data in annular bins and estimate the az- 
imuthally averaged VolU) from the visibility moment Vo(U) 
and VxiU) from Vi{U) and v^{U) in the same bins. 

• Assuming the power spectrum of the turbulent velocity 

fluctuations to follow power law, we consider ~ 

to model its power spectrum. We convolve 
with VoiU) to make the azimuthally averaged model 
'Pxmodi^)- important to realise here that as discussed 
in Section (3) our model of (U) derived this way does 

not have the proper amplitude. We shall restrict ourself in 
estimating the slope of the turbulent velocity power spectra 
here. The solid line in Figure (2) shows the column density 
power spectra for the case (b). The dotted line corresponds 
to VflU) for /3 = —2.0 while the dot-dash line shows the 
model "^ith 7 = —2.5. As expected, the model 

and the estimated value are in mismatch here. 

• We vary 7 and compare 'Pxmodi^l with Vf{U) by a 
least square estimator with likelihood function of the form 
exp(—L), where L is defined as 


L(7) 


LV) 

Lmin 


LV) = I] 


V 

' V 




(14) 


Here the index ‘i’ refer to a particular annular bin and N 
is the total number of bins. Note that both the power spec¬ 
tra are normalised at a given baseline Uo and hence this 
estimator is not sensitive to the amplitude of the velocity 
fluctuation spectra. We find the best fit value of 7 , i.e, j3 
and also error in estimating this value. 


Figure (3) summarises the efficacy of the model depen¬ 
dent maximum likelihood estimation of the slope of the tur¬ 
bulent velocity fluctuation power spectrum from the four 
sets of simulations. Let us concentrate on the panel (a) first. 
Here /3 is plotted against the different values of input jS for 
a = —1.5 and i = 10° with 3cr error bars. The solid line 
corresponds to an exact match between the input and esti¬ 
mated values of p. 

Clearly, for a lower inclination angle, i.e i = 10°, the 
estimator does a better job in accessing the real value of 
p. For the higher inclination angle, i.e i = 45°, for all the 
values of P the best estimate has systematically higher un- 
certaininity. This is expected as for high inclination angle 
our assumption of independency of Vz on line of sight coor¬ 


dinate is poorly valid. We also see for P values < —3.0 the 
measurement uncertainty is considerably high. 


5 DISCUSSIONS AND CONCLUSIONS 

In this letter we define visibility moments and an estima¬ 
tor of turbulent velocity fluctuations based on the visibility 
moments.Using numerical simulations we find that for lower 
inclination angles and turbulent velocity fluctuation power 
law slope > —3.0, the visibility moment estimator can be 
used. Though we have used a model dependent approach 
here, in principle it is possible to deconvolve the column 
density power spectra from Vf{U) and measure the veloc¬ 
ity power spectrum of the line of sight velocity fluctuations 
directly. Moreover, applying this method to the thin disk 
case, we can recover the turbulent velocity power spectrum 
completely. This would be particularly useful to access the 
energy involved in the large scale fluctuations seen in exter¬ 
nal galaxies. 

We have considered several simplifying assumptions. 
Firstly, we ignored the uncertainty in reconstructing the ro¬ 
tation velocity of the galaxy from the data, observational 
noise and incomplete baseline coverage. Each of these are 
expected to provide additional uncertainty for a real data 
and need to be looked upon carefully. Particularly, recover¬ 
ing the turbulent velocity power spectra depends highly on 
the accuracy of the rotational velocity estimate. For galaxies 
with warp disks, estimating rotational velocity is tricky and 
this method may not be useful there. Random fluctuations 
that may arise because of different systematic effects in the 
galactic rotation, however, are not expected to have scale 
invariant nature. An unwanted mixer of these fluctuations 
with the fluctuations generated because of the underlying 
turbulence would destroy the power law nature of the spec¬ 
trum. Hence, estimate of the velocity spectra itself can be 
used as a self diagnostic test. 

The approach described here to estimate the line of 
sight turbulent velocity fluctautions from the external galax¬ 
ies, apparently follow the same philoso phy from the method 
of st atistics of centroids of velocities (lEsauivel fc LazarianI 
l2005ll . where the unnormalised first moment of the spe¬ 
cific intensity function is used to estimate turbulent velocity 
statistics. Main difference in our approach is that we use 
the visibility function directly. Hence, it does not have any 
unwanted noise bias like that arise because of image recon¬ 
struction from incompletely sampled visibilities. 

Scaling relations of compr essible and highly magne ¬ 
tized turbulence are studi ed in iKowal fc LazarianI ll2007h . 
lEsauivel fc LazarianI 1 I 2 OO 5 II has also discussed incompatibil¬ 
ity of the velocity centroid method for high Mach number 
turbulence. Though the 10 kpc scales coherent structures 
seen at external galaxies are most likely be the result of 
ISM turbulence, it is neither clear if the neutral medium 
(H I ) at these scales would have an MHD coupling or they 
are highly supersonic. In either case, the applicability of the 
method described here needs to be investigated. Another 
basic assumption in developing this method has been that 
the H I is optically thin. The interpretation of the intensity 
spectra would cha nge in case the optical depth of the H i is 
sufficiently higher (iLazarian fc Pogosvanll2004h . 

Nevertheless, here we discussed an effective way of mea- 
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suring the statistics of the turbulent velocity fluctuations 
from external galaxies at large scales. We plan to use this 
estimator to the spiral galaxies with lower inclination angle 
from the THINGS sample in near future. 
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